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Q^' The interaction of classical gravitational waves (GW) with matter is studied 



^^. within a quantum mechanical framework. The classical equations of motion in the 



long wave-length limit is quantized and a Schroedinger equation for the interaction of 
\r2 \ GW with matter is proposed. Due to its quadrapole nature, the GW interacts with 

j^ ■ matter by producing squeezed quantum states. The resultant hamiltonian is quite 

different from one would expect from general principles, however. The interaction of 
GW with the free particle, the harmonic oscillator and the hydrogen atom is then 
studied using this hamiltonian. 
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§1. Introduction. 

When the two words "quantum" and "gravity" appear together in the same sentence, 
they usually appear as "quantum gravity" , a term which immediately brings to mind energies 
and length scales at the level of the Planck scale. With the advent of LIGO, and the prospect 
of the direct detection of gravitational waves (GW) [|T| , |]2[ , however, the interplay of classical 
and quantum gravity at a quantum mechanical level begins to take on added importance. 
Indeed, from the classical analysis the effect of GW on matter is expected to be so small 
that one begins to expect that a quantum mechanical treatment of this interaction is needed 
0, 0. It is, moreover, at the quantum mechanical level that experimental evidence for the 
interplay between gravity and quantum mechanics will first appear, rather than at the level 
of full quantum gravity. 

Surprisingly, the interaction of GW and matter at a quantum mechanical level has not 
been systematically explored. We know of two approaches currently in the literature. The 
first is by DeWitt [Q and is concerned solely with the general formulation of quantum 
mechanics on a curved, background space(time). The second is by Weber [^] and, as he 
was concerned with the response of GW detectors to incident GW's, deals exclusively with 
the interaction of GW with matter. His approach involved the direct quantization of the 
linearized classical equations of motion for a test particle interacting with a GW propagating 
on a flat background. Surprisingly, these two approaches are not equivalent to one another. 
One cannot, for example, obtain Weber's result, or the results of this paper, from DeWitt's 
general formalism by taking the linearized gravity limit. Fundamentally, we shall find that 
this is due to the different approaches taken by Weber and DeWitt. Weber begins with 
essentially the geodesic deviation equation which he then quantizes while DeWitt does a 
straightforward generalization of the quantization proceedure in fiat spacetime to curved 
spacetimes. 

In this paper we shall present a systematic study of the interaction of GW with matter 
at a quantum mechanical level in the long wave-length and small velocity limit. Like Weber, 
we shall start with the classical equations of motion, but unlike his approach we shall not 
first linearize the equations of motion. We find that the hamiltonian derived by Weber is 
valid only under very restrictive circumstances and, in particular, is not valid if the test 



particle is charged. In fact, due to the quadrapole nature of the GW (instead of the dipole 
nature of the electro-magnetic (em) field), the net effect of an incident GW on any quantum 
mechanical system is to produce a squeezed quantum state. This aspect of the interaction 
could not be seen within the linearized approach of Weber. Another aspect is that the GW 
is found to couple with test particle in a way that is quite similar to the minimally coupled 
em-field. While on the one hand this is what we would expect if the analogy between GW 
and em-waves is to hold, on the other hand it is quite different than what one expects from 
general principles. Spin-0 scalar particles like the ones we shall be concerned with in this 
paper are not expected to couple to the GW in this way. 

The rest of this paper is organized in the following manner. In §2 we shall quantize the 
classical theory using the usual canonical quantization proceedure. Then, in the succeeding 
sections §3 to §5 we shall use this hamiltonian to study the interaction of the GW with the 
free particle, the two dimensional harmonic oscillator, and the hydrogen atom. Concluding 
remarks and comparisons between DeWitt's, Weber's and our results can then be found in 

§6. 

§2. Quantum Mechanics. 

The equations of linearized gravity on a Minkowski background are well known. From 
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r^/3 = o^'^'^ (5a V + 5/3/^- - 5-^a/3) , (1) 



where h^^y is the perturbation off r/^j,, the flat Minkowski metric and 

^afiiiv = - {dfidphau + dadyhp^ — dad^hp^ — dpdyhan) ■ (2) 

We shall only be concerned with GW's propagating on a flat background and its effect on 
scalar, spin-0 particles. By choosing the transverse and traceless gauge for the GW, 

V = o, v^v = o, K = Q, (3) 

all the gauge freedom has been removed and the GW contains only its two physical polar- 
izations. The flrst polarization state hu = — /i22, ^33 = ^12 = ^21 = is usually called 
the +-polarization while the second hi2 = /i2i, hjj = (no sum) is usually called the 
X -polarization. In this gauge the only non-zero componants of the curvature tensor are 



^0,M- Q^ - 2 5t2 ■ ^^> 

As usual, greek indices shall run from — 4 while latin indices run from 1 — 3 and we are 
following the sign convention in 0. 

Classically, the response of a scalar test particle to the passage of a gravitational wave 
is given by the geodesic deviation equation. The overall effect of the GW is to introduce an 
effective tidal force to the equations of motion, 

^^ = -^Kmo^' + F' > (5) 

where x^ is the location of the test particle, m is its mass and Fj represents all other forces 
acting on the particle (see |^ for a complete derivation and explanation of the above). In the 
derivation of eq. (^, certain important approximations were made. First, it was assumed 
that the test particle is slowly moving and any velocity terms in the geodesic deviation 
equation may be neglected. Second, the long wave-length limit was taken, meaning that the 
reduced wave-length A/27r of the GW was assumed to be very much larger than the range 
of motion of the test particle. The GW can then be treated as a function of time only. 

Eq. d^), the classical equations of motion for the test particle, is our starting starting 
point. As we shall see, there are two different ways of quantizing the classical theory which 
are unitarily equivalent to one another only under very restrictive circumstances. We shall 
start with the most straightforward quantization procedure. For convenience, we shall as- 
sume that any external forces are time and velocity independent and can be represented as 
a potential V . Moreover, we shall always treat the GW as an external, classical field. 

The classical lagrangian for eq. (^ is 

X Tfl 

Lnr = -^mx] - —R^^qX^x'' -V , (6) 

where the dot denotes derivative with respect to time. Then choosing x^ as our generalized 
coordinate, the canonical momentum is pj = mx^ the mechanical momentum, and the 
hamiltonian is 

Hw = ^ + "^Ki^o^' + V. (7) 

We may then "quantize" this hamiltonian in the usual manner by replacing x^ and pj with 
operators which satisfy the canonical quantization condition: [x-^, p/c] = ih6l.. It is essentially 



a linearized version of eq. (|^ which is used by Weber and others [0 in studying the quantum 
mechanical properties of GW detectors. 

The second way to proceed is to note that the lagrangian eq. @ is equivalent to 

L = -mx^ — mTQf^XjX^ — V , (8) 

up to an integration by parts. If we once again choose x^ as our generalized coordinate, 
we now find that the canonical momentum is Pj = mxj — itlTqi^x^ while the hamitonian 
becomes 

H=^(p, + mTi,x')" + V. (9) 

Canonical quantization then gives: [x^ ,Pk\ = ih5\. The Schroedinger equation is then 

-ih^-^ = Hi; , (10) 

with the expectation value of any operator O defined as 

(V'|0|V)= fi^Oipd^x. (11) 



A priori, we might have expected the presence of Jg^^^ in the integral eq. (0) where (7*^^^ 
is the determinant of the three dimensional metric for a hypersurface. Since, however, we 
are dealing with linearized gravity, g^jj = 6jk + hjk and g^^^ ;^ 1 + h'jj^/2. To lowest order, 

Since Fg^ depends only on time, the quantization condition [x^,Pk] = iMl implies that 
[x^ , mxk] = ih6l also. The two hamiltonians eq. (0) and eq. (|) are related to one another by 
a simple canonical transformation. Thus there is a time dependent unitary transformation 
which maps one to the other. On this level one can just as well use either eq. (^ or eq. (|^) to 
analyze the response of the test particle to a GW. How the system is viewed physically, on 
the other hand, differs tremendously from one hamiltonian to the other. Notice the striking 
similarity between eq. (^ and the hamiltonian for a charged particle minimally coupled to 
the em field. Like the em-field, the GW is expected to carry momentum, which is recognized 
by eq. (^, and like the vector potential for the em-field, eq. (^ implies that the connection 
Fq^ plays a fundamental role at the quantum mechanical. One should keep in mind, however, 
that we are dealing with a scalar test particle and for a scalar particle one does not expect 



the connection to appear in the hamiltonian in this manner. Instead, at the most one would 
expect a direct couphng to the metric of the spacetime, as was obtained by DeWitt. 

Let us emphasize that this simple relationship between the two hamiltonians holds only 
under very restrictive conditions. For example, if the test particle is also charged, then under 
minimal coupling eq. (|^) becomes 



H= — [P,-^^A,+mTl,x'] +V, (12) 



Hw = ^{p,-lA^ +\%,,,x' + V, (13) 



in the non-relativistic limit while 

2 

2m V""^ c''V ' 2^ 

where q is the charge of the particle and Aj is the vector potential. Expanding out the 
quadratic momentum term, we find that H contains a direct interaction term between the 
em- field and the GW while Hy/ does not. Moreover, due to the presence of the em- field the 
canonical momentum Pj and the mechanical momentum mxj are no longer related to one 
another through a canonical transformation. The two hamiltonians will, in general, describe 
different physics. 

Actually, this type of equivalence up to a canonical transformation between two hamil- 
tonians also occurs when an em-field interacts with a charged particle. If an em-field in the 
long wave-length limit interacts with a slowly moving charged particle, then if we neglect 
any velocity terms in the interaction the hamiltonian can be approximated by 

He^ = ^ + qE,{t)x,{t). (14) 

where Ej is the electric field and is treated as a classical field. Because of the long wave- 
length approximation, is a function of time only. The charged particle actually couples 
minimally to the em-field, of course, but in this limit the two hamiltonians are related to 
one another by unitary transformation. 

We should also mention that both eqs. (|^) and (^) hold only in the long wave-length limit. 
Neither one is expected to hold in the general case where Fq^ is also position dependent. 
In fact, if we arbitrarily allow -Rgfco ^'-' ^^ position dependent, then the equations of motion 
calculated from eq. (|^) will contain terms which depend on the derivative of Rq ^q which 
are not present in the geodesic deviation equation. Indeed, we see that since the geodesic 
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deviation equation depends on the curvature, the hamiltonian calculated from it can only 
depend on the connection and should therefore have a form which is similar to eq. (^ and 
not to eq. (0). We thus conclude that eq. @ and its generalization eq. ([I2D are the correct 
hamiltonians for the interaction of a scalar test particle with a GW. 
§3. The Free Particle. 

In this section we shall study the interaction of a GW with a free test particle using the 
hamiltonian in eq. (^. The solution of this problem is fairly standard (see, for example ^). 
In general terms, due to the quadrapole nature of the GW the effect of an incident GW is 
to produce a squeezed quantum state, which is well known from quantum optics (see for 
a complete review). Nevertheless, we shall present a fairly complete analysis of this system 
as the analysis of the interaction of the GW with the two-dimensional harmonic oscillator 
follows in much the same way. 

We first choose the z-axis to lie parallel to k, the wave-vector for the GW. Then, due to 
the transversality condition eq. (^), Fq^^ has non-zero componants only in the x — y plane 
and the particle undergoes free motion along the z-direction. The problem reduces to that 
of two dimensional motion with the hamiltonian 



p2 

H = ^+ nkX'Pk + ^rikni^'^' , (15) 



where we have used the traceless condition and the indices j and k now run from 1 — 2 only. 
Because we are dealing with linearized gravity, the last term in eq. (|T3]) dependent on F^ 
has little meaning and we shall set it to zero by hand. 

Oftentimes we shall find it more convenient to work with 2x2 matrices representing the 
two possible polarization states of the GW. When doing so, the GW is witten as a matrix h 
with the explicit componants of h being hjk for j,k = 1,2. Moreover, using the symmetry 
properties of the GW we can express 

h,k = 2/(t) (e^a], + e+cT%) = 2f{t)eAaf, , (16) 



where a^ and a^ are the first and third Pauli matrices, respectively and the index A runs 
from 1 — 3. 2/(t) is the amplitude of the GW while ^x and e+ are the polarization states of 
the wave. They are in general time dependent and satisfy 

4+4 = 1, (17) 



for all t. We shall also assume that the GW incidents the particle at t = so that /(t) = 
for t < 0. 

We next define raising and lowering operators 

^' = (^) («. + 4) ' P3 = -^ (^) («. - 4) , (18) 

where w is, as we shall see, related to the uncertainty in the initial position and momentum 
of the particle. Then, 



H 
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[ayaj + 2 j ~ X V ^' ^ "■'■ j " J^iH^) (ajafe - SO/lj • (19) 



Working in the Heisenberg representation, we find that 



daj w i 1 

dt 2 ^ ^ ^^2 



-77 = -^ — [aj - a-j) + -hjkak ■ (20) 

Next, note that the raising and lowering operators must satisfy the commutation relations 

[a,(t), ak{t)] = , [a jit), 4(t)] = 5jk , (21) 

at equal times. This implies that aj{t) are related to aj(0), the free operators at t = 0, by 
the canonical transformation 

«i(^) = UjkakiO) + Vjkal{0) , 

a]{t) = al{0)ukj + ak{0)vkj , (22) 

with the bar denoting the complex conjugate. Eq. ( P^ ) is a time-dependent Bogoluibov 
transformation with u and v being the generalized Bogoluibov coefficients. They are 2x2 
complex matrices which, due to eq. (pT]), must satisfy 

Mf * = u^v , mm' — vv^ = I (23) 

written in matrix form, t denotes transpose and / is the identity matrix. Moreover, they 
have the boundary conditions u(0) = I and f (0) = 0. Then, from eq.(^), 

— -(m + f^) = —izuiu — v"^) H — (u + v^) , 
dt 2 

^^(u - v^) = ~{u - v^) . (24) 



Eq. (^) is difficult to solve for general hjk- For the special cases of linearly and circularly 
polarized GW, on the other hand, we can find exact closed form solutions to eq. (P^). In the 
case of linearly polarized GW, the polarization states e^ are constant while f{t) is arbitrary. 
One can then directly intergrate eq. (|23) to obtain 



ui = I cosh/(t) -2y[ cosh/(t) f\osh{2 f{t'))dt' 

-sinh/(t) l\mh{2f{t'))dA \l 

-i-|sinh/(t) f\osh{2 f{t'))dt' 
2 [ JO 

-cosh/(t) l\mh{2f{t'))dt'\eA(r^ , 
vi = z^l cosh/(t) j\osh{2 f{t'))dt' 
-sinh/(t) rsinh(2/(t'))t^t'|/ + 

1 sinh f{t) +ij( sinh f{t) t cosh{2 f{t'))dt' 

- cosh f{t) t sinh(2/(t'))rft' ) [ca^^ • (25) 

For circularly polarized GW, on the other hand, f{t) = /o, a constant, for t > while the 
polarization states now varies with time 

^ - "- • ^ - -"'- ■ (2«' 

The exact solution to eq. (|^) in this case is quite complicated and we shall present here 
only the small /o solution 

Ucir = \^ 7) 1" "iT" (cosh ilt — COS fit) > / 

< cosh fit + COS fit — 2 — iwt + i — sin fit > e • a 

2 I fi ^ 

H {sinh fit — sin fit} — — — 

+— {sinh fit + sin fit - 2fit} i^^r^ ■ ^ , 
y^.^ = J '^— + l^ (cosh fit + cos fit - 2) W 



< cosh fit — cos fit — 2 + itx^t — i-- sin fit > £* ■ a 

2 1 fi J 

H {sinh fit + sin fit} — — — 

-^{sinh fit-sin fit} i^^ -a, (27) 

where we have used the vector notation defined in Appendix A. We refer the reader to 
Appendix A for the derivation of the exact solution. 

To complete the quantum mechanical analysis, we note that the unitary evolution of 
aj{t) can be implimented by the following unitary transformation on aj{0) and aKO), 



aj, 



,j(t) = Rit)S{t)aj{0)S\t)R\t) , a]{t) = R{t)S{t)a]{Q)S\t)R\t) , (28) 

where 

S{t) = exp{(p,fea,(0)a,,(0) - p,fca](0)4(0))/2} , (29) 

is the generalization of the squeeze operator of quantum optics, while 

R{t) = exp{-tejka]{0)ak{0)} , (30) 

is the generalized rotation operator, p and e are 2x2 matrices with p symmetric: p = p*. 
It is then straightforward to varify that S is indeed unitary while R is unitary if and only if 
e = e"!" is hermitian. The time evolution operator is then U{t) = R{t)S{t). 

Using the canionical commutation relations, we find that u and v must be related to p 
and e through 



u={J2 "-^^^^ 1 e" 



-fET^Ve-. (31) 

The constraints eq. ( pSj) requires that e also be a real matrix, while [p, p^] = 0. p and e are 
the generalization of the squeeze parameter and the rotation angle of quantum optics 
with p containing within it the squeeze angle also. 

The system has now essentially be solved and one need only specify the initial state 
of the particle. This involves specifying not only the initial average position {x^{0)) and 
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momentum (p'iO)) of the particle, but also its initial uncertainty in either position and 



momentum. Since Ax(0) = Jh/{2mw) and Ap(0) = Jhm-cu/2, doing so determines w. 
§4. The Harmonic Oscillator. 

We now consider a GW incident on a two dimensional harmonic oscillator. Once again 
for simplicity we shall assume that GW travels in a direction perpendicular to the oscillating 
plane of the oscillator, so that 

H = ^ + \mu'x] + Ti,x'p, , (32) 

where uj is the oscillation frequency of the pendulum, m is its mass and j. A; = 1, 2. The oscil- 
lator is assumed to have no equilibrium length in distinct contrast to the system considered 
by Weber. We have once again neglected the F^ term. 

We again define raising and lowering operators as in eq. (0), but with uj instead of vj, 
and obtain 

H = fiuj{a},aj + 1) hjk [cLjak — alal) • (33) 

Note that when e^ = 1 and 6^=0, the x and y directions decouple and we have the 
hamiltonian for a one mode squeeze state (see 0]) in each direction seperately. When, on 
the other hand, £+ = and £x = 1, the two directions couple and we have the hamiltonian 
for a two mode squeeze state. 

Proceeding as before, we look for solutions of the Heisenberg equations of motion with 
the form eq. (|2^). We again obtain a set of matrix differential equations for u and v, but 



unlike the case of the free test particle, these equations are no-longer easily solveable even 
in the case of linearly and circularly polarized GW. 

Let us first consider the case of linearly polarized GW for which hjk = 2f{t)eAO'fk = 
2f{t)ajk and Ea a constant. The equations we now have to solve are 



—- = —iwu + fav^ , — - = —iwu + jaw , (34) 

dt dt ^ ^ 

along with their corresponding complex conjugate equations. Although under certain in- 
stances these equations can be solved exactly in terms of Matthieu functions, doing so will 
not be physically illuminating. Rather, we note that usually \f{t)\ -C 1 and we can solve 
eq. (|34D perturbatively about its /(t) = solution 
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e'^^'ll + L^ + Aco^ P f{t') t' f{t")e^'^^''-'"Mt'dt" 
\ 2 Jo Jo 

+2iuof{t) t /(t')e'*"(*-*')rft - 2iuj r f{t')dt'\l, 

f{t) - 2%uj r /(t')e"''"(*-*')rft'') ~G . (35) 



As we have mentioned, the hamiltonian for the interaction of a GW with a harmonic 
oscillator is very similar to that of the hamiltonian for squeezed states in quantum optics. 
Indeed, if we solve eq. ( PID in this case, we find that p = re^^'^a and e = el where 

tan(e - wt) = 2ujf{t) [ f{t') cos(2cu(t - t'))dt' -2uj ( f{t')dt' 

Jo Jo 

+4uj f' fit') t f{t")sm{2uj{t'-t"))dt'dt" , 
Jo Jo 

, ( 0^^ A 2u;J^f{t')sm{2u;{t-t'))dt' 

ta.n(e — 26 + wt) = ^-- — -7 — — — . (36 

^ ' fit) -2u; JU if) sm{2uj{t-t'))df ^ ^ 

There is, however, some ambiguity in the solution for r. If one uses the equation for u, one 
finds that 

r' = fit) - iufit) f fit') sin(2c.(t - t'))dt' 
Jo 

rt ft' 

, fit') 

Jo Jo 

while if we use the equation for v 



+8cu2 /* fit') t fit') cos(2^(t' - t"))dt'dt" , (37) 

Jo Jo 



2 



r2 = /2(t) _ 4ujfit) I fit') sin(2cu(t - t'))dt' + 800^ ( f fit') cos(2w(t - t'))dt'\ . (38) 

This ambuigity araises from our perturbation scheme and the desire to polar decompose u 
and V into a modulus and phase. The differences between the two expressions only occur 
in the last term, however, and are very small for wt <^ 1. Any ambiguity disappears if 
we use eq. (^) in eq. ( ^4|) to obtain differential equations for p and e and then preform 
the perturbative analysis on p and e directly. For our purposes, however, the perturbative 
solutions eqs. (^) and (|37D are sufficient. 

The parameters r, e and (p are usually called the squeeze parameter, the rotation angle 
and the squeeze angle, respectively. There are two properties of r and (p that are of interest. 
First, from eq. (|35|) we expect a resonance to occur when the frequency of the GW is equal 
to 2uj. This can be seen explicitly by taking fit) = fosinQt for a monochromatic GW 
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and performing the integral in eq. (|35|) explicitly. Physically, this is due to the quadrapole 
nature of the GW, and its quadratic coupling to the harmonic oscillator. At this order of 
the perturbative expansion we get frequency doubling, and thus a resonance at i7 = 2u; 
for u and v. If we were to then go to the next perturbative order, we would get frequency 
quadrupoling and a resonance at Q = 4uj and so on. 

Second, notice that to this order (p is independent of the amplitude of the GW. Conse- 
quently, a very weak GW can still produce a large response in the squeeze angle. Indeed, 
for the monochromatic wave f{t) = /o sin fit, 

u sin Qt — Q sin 2ujt , , 

tan(e - 20 + wt) = -2uj-— -, , (39) 

^ ^ ^ (fi2_2cu2)sinfit-2cufisin2cut ' ^ ^ 

when Q ^ 2uj. Note that the denominator of eq. (^) may vanish for certain t. (This 
will, in fact, almost certainly happen when if the incident GW is on resonance Q = 2uj.) 
Consequently, we can expect very large variations in the squeeze angle no matter how 
weak the incident GW is. Since, however, the squeeze parameter r ~ f{t), this large phase 
shift will still be difficult to see. Note also that this does not happen for e, which depends 
on /2 and to lowest order is unaffected by the GW. 

Finally, let us consider a circularly polarized GW for which 

dU . r ■ A i ^^ ■ r ■ A i f Ar^\ 

— - = -twu + foe ACT f ' , -J- = -iwu + foSAO- m' . (40) 

at at 

Although these equations can be solved exactly using the techniques described in Appendix 
A, the final result will not useful as it involves the solution of a fourth order polynomial 
equation. We shall instead once again solve eq. ( ^OD perturbatively and obtain 

Uair = e -'"*<! 1 - 2iijfo t e{t) ■ e{t')e^'^^'-''^dt' 

+4u;2 lY' e{t')e{t")e^'''^''-'"^dt'dt"\l 

il Jo ii li 

.2z^ [''Eit)-s-{t')e^'-^'-''^dt' 

il JO 




+^^^, X Jo ( "^^^ ■ "^^'^"^^^ ■ "^^"^ 
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e ■ a 



Voir = /oe'"* (l + 2iLU I e{t)e{t')e-^'^^'-''^dA 

+2^/0^6^'^* (/o* *) ■ ^it'y-''''^'-'^dt) ^ • (41) 

We once again expect a resonance at fi = 2uj. The corresponding p and e can also be found 

in this case. 

§5. The Hydrogen Atom. 

As we shaU be using time dependent perturbation theory to analyze the interaction of 
GW with the hydrogen atom, in this section we shall keep to generalities. An exposition of 
time dependent perturbation theory can be found in any quantum mechanics textbook (see 
for example ||lOl)- From eq. (0) the relevant hamiltonian is 






--Ti,A,x' + Fo.^V, + i^n^Tl^x^x^ (42) 

where p is the reduced mass of the hydrogen atom, Aj is the vector potential and V is 
the usual Coulomb potential for the hydrogen atom. Although the vector potential will be 
treated as a field operator, the energy density for the em-field has not been included in 



eq. (0). We have moreover neglected the spin of the electron, as we have not yet considered 
the effect of the GW on the spin of a particle. We also note that even at atomic energies, 
the wave-length of the GW will still be much larger than the size of the hydrogen atom and 
the long wave-length approximation is still valid. 

The first four terms in eq. (^2|) are what one usually obtains for a minimally coupled 
hydrogen atom. The additional three terms comes from the interaction of the GW with the 
atom and we shall treat these terms perturbatively. As usual, since the last term involves 
r^ we shall not consider its affects. 

Consider first the interaction term 

H^t = ^hjkX^ , (43) 

which couples the GW directly to the electron and can cause either the absorption or emis- 
sion of a "graviton" . From time-dependent perturbation theory the transition probability 
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between an initial state \init) and a final state \fin) is 

2 2 

|c(+oo)p = -—i^ujtf\hjk{oJif)f'\{fin\ (xjXk - 5jkX^ /i) \init)\^ (44) 

where hjk{wif) is the fourier transform of hjk{t) and huif is the energy difference between the 
initial and final states. The quadrapole nature of the interaction can now be seen explicitly. 
From the Wigner-Eckhart theorem the allowed transitions are: Am = ±2, A/ = ±2 where 
m and / are the angular momentum quantum numbers. 

Using eq. (||) and the fact that the intensity of the GW is / = /i^fc)u;^c^/327rG, find that 
the cross-section for this interaction has the form 

<ycross{^) = Jifll , (45) 

where J,/ is a numerical factor depending on the initial and final states and Ip = Gh/c' is 
the Planck length. (The reader should not confuse this cr cross with the Pauli matrices.) The 
cross-section is thus proportional to /^ = 10~^^ cm^, and is extremely small no matter what 
value Jif takes. Combined with the expectation that the intensity of GW's in the universe 
which are at atomic energies are very low, such a small cross-section means that it will 
essentially be impossible to see a GW induced transition in a hydrogen atom. Nevertheless, 
it is interesting to see the appearance of the Planck scale even at this simplestic level. Notice 
also the absence of knowledge of the em nature of the system, such as the charge of either 
the electron or the nucleus, in eq. (^). This once again underscores the fact that HlJ^ is 
due to the direct interaction of gravity with matter. 
Let us next consider the interaction term 

H^, = ^h,uA,x^ . (46) 

In may ways this is much more interesting as it couples the GW with both the em-field and 
matter. In particular, this is a dipole coupling in contrast to the quadrapole coupling of 
ifjj. The absorption of a graviton can thus cause the emission of a photon along with the 
excitation of the hydrogen atom to an excited state. Conversely, the emission of a photon can 
also cause the emission of a graviton. Theoretically, the passage of a GW should shorten the 
lifetime of the any excited state of the atom. Indeed, if we combined eq. (^) with the usual 
dipole interaction term in eq. (|4^) , we find that the cross-section for emission or absorption 
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of a photon in the dipole approximation will be slightly larger by a factor 1 + const. h'jf,{0) 
where the constant depends on the polarization state of the GW. Since, however, hjk is 
extremely small, this shortening cannot be seen experimentally. 

Using Fermi's Golden Rule, we can also calculate the transition cross-section for the 

(2) 

emission or absorption of a single photon due to HIJ^ , 

(^crossii^) = J'ifOi {uJifabohr/cf fp , (47) 

where once again J[f is a numerical factor depending on the initial and final states of the 
hydrogen atom, huif is the energy of the photon, ai^ohr is the Bohr radius and a is the fine 
structure constant. For typical atomic transitions, (Tcross ~ 10"^/^ which is even smaller 
than that of the direct interaction Hl,J^ since the transition is mediated by the emission of 
a photon. 
§7. Concluding Remarks. 

To conclude, we have quantized the classical equations of motion for a scalar, spin-0 
interacting with a classical GW in the long wavelength limit when the velocities of the 
particle is small. We find that due to the quadrapole nature of the interaction the effect of 
the GW is to produce a squeezed quantum state El. We have then used this formalism to 
calculate the effects of the GW on the free particle, the harmonic oscillator and the hydrogen 
atom. 

The hamiltonian H that we have obtained is quite peculiar, however. On the one hand 
it has the form one would expect if the analogy between the GW and the em-wave is to 
hold at this level. Like the em- wave, the GW carries momentum and energy. And like 
the em-wave, since the classical equations of motion depends on a physical observable (the 
field strength for the em- wave, the curvature tensor for the GW), the hamiltonian must 
depend on their respective connections. A coupling of the form {Pj + ttiFq^x'^)^ would thus 
seem quite natural. Moreover, like the vector potential for the em-wave it means that the 
connection Fq^ now takes on an independent meaning at the quantum mechanical level. On 



^Curiously, the creation of primordial GW from the de Sitter vacuum is also due to quantum 



mechanical squeezing [11 
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the other hand, this couphng resembles minimal coupling for the gravitational field. As we 
are dealing with a scalar particle, from general principles alone we would not have expected 
a coupling of this form. We would instead have expected a hamiltonian of the form given by 
DeWitt [§]. One cannot, however, obtain H directly from DeWitt's hamiltonian by taking 
the linearized gravity limit. 

DeWitt in |^ considered the quantum mechanics of a particle constrainted to move on 
a curved space(time). Then if gij is the metric on the three dimensional hypersurface he 
obtained a hamiltonian of the form 

for a non-relativistic particle where g is the determinant of gjk- a Let us now try to apply 
this hamiltonian to the case of GW. In this case gjk = 6jk + hjk and due to the traceless 
condition on hj^, (7 = 1 up to first order in hj^. Consequently, 

This hamiltonian is very different from only the one we have derived, but also from the 
one used by Weber. Moreover, from the transformation law for general coordinate transfor- 
mations (eq. (5.42) of 0), -ff_DVK cannot be mapped to H through a judicious coordinate 
transformation. 

We should also mention that in DeWitt's derivation of How he actually obtained pj = 
—ih{dj + dj ln5f/4) — dj^ where $ is an arbitrary real function oi Xj and t (eq. (5.24) of ||^). 
$ was then set to zero by arguing that it can always be removed by a local unitary phase 
transformation on the wavefunction. One finds precisely a term of this form in H. This, 
however, still does not explain the presence of the second term in How which is absent in 
H. 

We believe that the differences between How and H to be more fundamental than a 
simple gauge transformation of the wavefunction, however. DeWitt's derivation of How 
was based on a natural generalization of quantum mechanics on flat spacetime to curved 



^DeWitt considered motion on an n-dimensional space. We need not be so general. 
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spacetimes. Our hamiltonian, on the other hand, is based on the classical geodesic deviation 
equation. Physically, this is because the passage of the GW will affect not only the test 
particle, but also the observer (measuring apparatus, laboratory setup). Since the observer 
is also part of the universe, it cannot be isolated from any gravitational effects and also 
responses to the GW. Consequently, one cannot measure (observe) the absolute motion of 
the test particle, but rather the relative motion between the observer and the test particle. 
This motion is governed by the geodesic deviation equation. In DeWitt's formalism it is 
not certain where the observer is. Since he deals with only a single particle, we suspect 
that the observer has tacitly been "removed" from the spacetime, something one can do in 
Newtonian mechanics, but cannot do in general relativity. As the results of this paper is 
only valid within the long wave-length limit and for a fixed choice of gauge we have not yet 
been able to show this explicitly. 

Further study of the interaction of GW with matter is clearly needed to clarify these and 
other issues. In particular, a version of the hamiltonian which is valid outside of the long 
wave-length approximation is needed and the effect of GW on particles with spin should also 
be incorporated [|12| . Once this is accomplished, a thorough comparison of our formalism 
and DeWitt's results can be made. 

We end this paper by briefly mentioning the differences between the hamiltonian used 
by Weber and the one we derived. As the work by Weber was mostly concerned with 
the properties of the interaction of GW with gravitational wave detectors, he had at their 
disposal a macroscopic length /. For the beam detector this is the length of one arm of the 
interferometer, while for the bar detector this is the equilibrium length of the spring for a 
harmonic oscillator. In both cases, the hamiltonian that was used has the form 

2 
D Tfl 

Hweb = \- TT^^^^ + "^-^0 10^^^ ) (50) 

2m 2 
along, say, the x-direction. (The harmonic oscillator term is absent of the beam detector). 
Importantly, x is measured from its equilibrium position. If we replace x ^ x + I in eq. (^ 
and expand out the quadratic term, we obtain Hweb up to a term proportional to x^. This 
was dropped in comparison to xl since the response to the GW is expected to be quite small 
classically. While this is perfectly valid in classical dynamics, since x is now an operator 
doing so is somewhat questionable at a quantum mechanical level. We have, however, re- 



preformed the analysis of both the beam and the bar detectors using the full hamiltonian 
eq. (^ in a quantum mechanics framework and found very little differences in the results 
obtained from using Hweb verses H. Consequently, as long as the linearization of eq. (H) is 
valid, and as long as the test particle is not charged, the hamiltonian used by Weber and 
others gives a good description of the interaction of a GW with matter. 

APPENDIX A 

We first define ^ = u + v"^ and x = ^^ — ''^^ so that eq. ( p^ ) becomes 

_ = -itux + fQeA(y ^ , -^ = -foeACT x , (M) 

for circularly polarized GW. Next, we note that any 2x2 complex matrix M can be written 
has M = 9qI + 9Aa^ where 9q and 9 a are complex numbers. We next consider 9 = {9i,92, 9^) 
as being a vector in a three dimensional complex space. Clearly the polarization states of the 
GW can also be represented as a vector sin this space. Moreover, e, e and £*x eaxe mutually 
orthogonal and thus form a natural coordinate system for this space. Consequently, we look 
for solutions of ( [Al| ) with the form 

X = AI + Be-a + C^^ + Dt^ ■ a , (A2) 

where A, B, C, D can be complex functions. Then 

HA 

— + hnc = , 

at 
dB 

— -nc- fo^D = , 

at 

dC 

-— + nB + foQA = , 
dt 

^-f,nB = 0. (A3) 

Next, using the boundary condition x(0) = n(0) = /, we find that 
u"^ cos z/_a; — z/^ cos zz+x + /q (cos z/_a; — cos u^x) 



A{x) 



z/i — z/?. 



B[x) = — K 5 (cosz/_a; — cosz^+x) 

/o n - ^~ 
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C{x) 



1 z^_|_z^_(z/+ sin i'_x — z/_ sin i^+x) + fo{i^~ sin z/_x — z/+ sin i^+x) 

To ^i^^Z 






where x = Qt and 



1 



4 = 2 l^ " ^^' ^ V 1 - 4/o^ j • (A5) 

The solution for ^ follows in exactly the same way. We once again take 

^ = EI + Fe-a + G^^ + m^^-a, (A6) 

and now find that 

K^ cosh K_x — K^ cosh K+x + /q (cosh K+x — cosh fi;_x) izuC iwD 

F(x) = | K-/o)(^-^-/o) (eoshK+x - cosh/.„x) + '"^^ 



G(x) 



1 /£+/!;_(«;+ sinh k^x — k_ sinh k^x) + fo{K+ sinh k+x — k_ sinh k_x) 



/o ni- K 



2 ,^2 



^ ^ «-/o)(/t^-/o) / sinh/t+x _ sinh/t_x \ 



1 



with 

4^-(^l + 2/^±Vl+4/o^j. (A8) 

Notice that for /o ^ 2, both v± and fi;± are real. Eq. (^) now follows if we take the /o ^ 
limit and keep terms up to 0(/q). 
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